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A general algebraic procedure for constructing coherent states of a wide class of exactly 
solvable potentials e.g., Morse and Poschl- Teller, is given. The method, a priori, is poten- 
tial independent and connects with earlier developed ones, including the oscillator based 
approaches for coherent states and their generalizations. This approach can be straight- 
forwardly extended to construct more general coherent states for the quantum mechanical 
potential problems, like the nonlinear coherent states for the oscillators. The time evolution 
properties of some of these coherent states, show revival and fractional revival, as manifested 
in the autocorrelation functions, as well as, in the quantum carpet structures. 
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INTRODUCTION 



Coherent states (CS) and their generalizations, in the context of harmonic oscillators, are well- 
studied in the literature [1-4]. Algebraic approaches have been particularly useful for providing a 
unified treatment of these states and their inter-relationships. For example, in Ref. [5], not only 
£T) . a general procedure for constructing a large class of oscillator based CS has been provided, but it 

has also been shown that, some of these states are dual to each other in a well-defined manner. 
The algebraic approaches straightforwardly lead to the construction of squeezed and other states 
showing interesting non-classical features. These elegant and powerful algebraic procedures of 
construction owe their origin, partly, to the simplicity of the Heisenberg-Weyl algebra; [a, at] = 
^ 1, characterizing the harmonic oscillator. Based on the symmetries and keeping in mind the 

desired requirements, various procedures have been developed for constructing CS for Morse [6-9], 
Hydrogen atom [10], Poschl- Teller [11-14] and other potentials [15]. The role of Morse potential 
in molecular physics is well-known [16]. The study of CS for hydrogenic atoms have assumed 



increasing importance in light of their relevance to Rydberg states [17], which may find potential 
application for quantum information processing [18]. 

It is known that, all the criteria desired of a coherent state and found in the oscillator based CS, 
^ ■ e.g., minimum uncertainty product, eigenstate of the annihilation operator (AO) and displacement 

operator states are not simultaneously achievable in other potential based CS. Hence, a number of 
CS, diluting one or more of the above criteria, have been constructed in the literature. 

In the context of algebraic approaches, super symmetric (SUSY) quantum mechanics [19] based 
raising and lowering operators have found significant application. In particular, eigenstates of 
the lowering operator for Morse [6] and Poschl- Teller [14] potentials have been found and their 
properties studied. Recently Antoine et. al. [13] have constructed Klauder type CS for the 
Poschl- Teller potential, using a matrix realization of ladder operators, their motivation being the 
temporal stability of the CS. It is well-known that the SUSY ladder operators act on the Hilbert 
space of different Hamiltonians except for the case of the harmonic oscillator. Establishing a precise 
connection between the complete set of states, describing the above CS, and the symmetries of these 
potential problems has faced difficulties [8, 20]. To be specific, in case of the Barut-Girardello CS 
for the Morse potential, the ladder operators are taken to be functions of quantum numbers, which 
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has led to problems in denning a proper algebraic structure. To resolve the same, some authors [8] 
have resorted to the introduction of additional angular coordinates [21]. 

In this work, we provide an algebraic construction of the CS for a wide class of potentials, 
belonging to the confluent hyper geometric (CHG) and hypergeometric (HG) classes. The procedure 
is based on a simple method of solving linear differential equations (DEs) [22], which enables one 
to express the solutions in terms of monomials. In the space of monomials, it is straightforward to 
identify various types of ladder operators, their underlying algebraic structures [23], and construct 
lowering operator eigenstate in a transparent manner. The fact that, the monomials and the 
quantum mechanical eigenfunctions are connected through similarity transformations, enables one 
to preserve these algebraic structures at the level of the wavefunctions and simultaneously obtain 
the desired CS. Thus the coherent state is initially potential independent. The information about 
a specific potential is then incorporated by fixing the parameters of the series and also the ground 
state wavefunction of the potential under study. The known results for CS are obtained in specific 
limits. In addition, our procedure demonstrates the construction of more general CS, similar to 
the nonlinear CS [24] in the oscillator example. The origin of confusion in identifying the algebraic 
structure in SUSY based approaches is then pointed out and subsequently resolved in a natural 
manner. It is shown that the recently found CS for various potentials are related to the two 
different realizations of SU(1,1) algebra. 

The paper is organized as follows. Keeping in mind the connection of the wavefunctions of the 
exactly solvable potentials, to be considered here, and the solutions of the CHG and HG DEs, we 
first briefly outline the relationship between the solutions of the above equations and the space of 
monomials. This connection is then made use of, to identify the ladder operators and the symmetry 
algebras in the space of monomials. AOCS are constructed, through a recently proposed procedure 
for solving linear DEs. This method yields the precise connection between the present, and the 
earlier oscillator based approaches for constructing CS. In Sec. Ill, the potential independent 
CS, thus constructed, are connected with Morse and Poschl- Teller potentials, as examples. The 
advantage of constructing the ladder operators in the space of monomials is then pointed out 
by resolving the difficulties in identifying the symmetry generators in the SUSY based approach. 
Various earlier known CS are then derived as special cases. Time evolution properties of some of 
these CS are then studied. Autocorrelation functions and the underlying quantum carpet structures 
[25] of these CS clearly reveal the phenomena of revival and fractional revival [26]. We conclude 
in Sec. IV after pointing out a number of open interesting problems, where the present procedure 
can be profitably employed. Keeping in mind the rich symmetry structure of the hydrogen atom 
problem and the various procedures employed for constructing the corresponding CS, we desist 
from the analysis of these CS here; this will be taken up in a future project. 

II. ALGEBRAIC CONSTRUCTION OF COHERENT STATES 

This section is devoted to the construction of CS in a manner which is potential independent. 
We make essential use of the fact that the exactly solvable potentials, to be studied here, belong to 
the CHG and the HG classes, whose solutions can be connected to the space of monomials through 
similarity transformations. In the space of monomials, the identification of symmetry algebras and 
their ladder operators become easy. The AO eigenstates in the space of monomials are first found 
through the above mentioned procedure of solving DEs and then connected to those at the level 
of the polynomials, through similarity transformation. Connection with earlier oscillator based 
approaches is also exhibited. 

A single variable linear DE, of arbitrary order, can be cast in the form [F(D)+P(x, d/dx)\y{x) = 
0; where F(D) is a function of the Euler operator D = xd/dx and P(x, d/dx) contains rest of the 
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operators. A formal series solution of the above DE is given by [22, 23] 

1 

_ ir _ 

n=0 



F(D) 



P(x,d/dx) 



(1) 



provided the condition F(D)x a = 0, is satisfied. This condition fixes the value of a. Here C a is 
the normalization constant. The proof is straightforward and follows by direct substitution. In the 
CHG and the HG cases, a = n (integer) corresponds to polynomial solutions and for a / n series 
solutions are obtained. 

Following an earlier work [23], the simplest set of operators, at the level of monomials, which 
give rise to a SU(1,1) algebra can be written as [27] 



K + =x,K- = 



' x d 2 b J_ 

dx 2 dx 



, and_K"3 = 



d_ b 

X dx' + 2 



(2) 



Here b is deliberately taken same as the CHGDE parameter, so that the algebra is well defined at 
the level of the wavefunctions. It should be noted that, these symmetry algebras are not unique, 
and the interesting consequences arising from this fact will be pointed out later. The AOCS 
corresponding to K- satisfies, 



K-<p(x,0) = -P<p(x,P). 



Left multiplying by x: 



d 



d 



x z —, + bx— + fix ip(x, 0)=O 



dx 2 



dx 



(3) 



(4) 



one can identify F(D) = (D + b — 1)D and P(x,d/dx) = 0x. The condition F(D)x a = gives 
a = 0, 1 — b. For a = one obtains, 



n=0 

oo 



(D + b-l)D 
(-l) n {f3x) n 



x 



n=0 



n! 



(b)n 



exp 



(D + b-1) 



x 



(5) 



N(P) _1 is the normalization constant. The commutation relation, [D,x d ] = dx d , and the definition 
of the Pochhammer's symbol, (b + n — 1) • • • b = T(b + n)/T(b) = (b) n , have been used in the 
above derivation. We can make immediate contact with the oscillator based approach of Ref. [5]. 
Denoting 



(D + P-l) 



x, 



(6) 



it is found that, they obey the Heisengerg-Weyl algebra, = 1. Such conjugate operators 

were explicitly constructed in [5] and used to obtain a large class of multi-photon coherent states. 
Knowing that the solution of the CHGDE can be written in the form [22] 



<&(-n; b; x) = (-!)" 



T(b) 



T(b + n 



-e~ K -x n . 



(7) 



one can obtain the coherent state at the level of the polynomial and later on at the level of 
the wavefunctions through appropriate similarity transformations, as illustrated below. Denoting 
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S = exp(— KJ), one obtains SK^S 1 S(p(x,f3) = (3S<p(x,(3) = /3<p(x,(3), where <p(x,(3) is the 
coherent state at the level of the polynomial: 

<p(x, f3) = NiP)- 1 ]T b; x) . (8) 

n=0 

It will be seen explicitly in the subsequent section that introducing appropriate measures through 
similarity transformations one would obtain CS for various quantum mechanical systems. It is 
worth emphasizing that the algebraic structures are transparently preserved in this approach. 
Analogously for HG case, taking the lowering operator to be [23] 

1 ( d 2 d \ 

K - = TDTb)[ x d^ + c Tx)> (9) 

the coherent state at the level of polynomials is given by 

X(x )7 ) =N( 7 )- 1 Y / 1 j2F 1 (-n,b;c;x), (10) 

n=0 



since 



We can connect to the oscillator based approach in this case also; calling 

+ (D + c-1) V ; 

it satisfies the oscillator algebra: [K-,K + ] = 1. It should be noted that the states obtained using 
the lowering operator K_ are nonlinear CS [24]. They are defined to be eigenstates of the AO of 
the type /(jV)T_, f(Af) being a function of the number operator in the oscillator case and the 
Euler operator in the present one; T_ is an arbitrary lowering operator. It was recently shown [28] 
that, the AOCS and the Perelomov CS are unified in the framework of nonlinear CS. 

Equations (8) and (10) reveal that these CS are quite general; since the eigenfunctions of the 
exactly solvable potentials to be considered here, arise as special cases of CHG and HG series, their 
corresponding CS will also follow from the above two general expressions. In fact, starting from 
the lowering operator 

( D + bl )... {D + bp ) d 

{D + ai) ■ ■ ■ (D + a p ) dx ' K ' 

more general nonlinear CS can be found: 

£(x, a) = 22 — p+ iFp(ai, • • • , a p , -n; b u ■ ■ ■ , b p ; x). (14) 

UTh. 

Here p +\F p is the generalized HG series. For the sake of completeness it should be pointed out, 
the above summation yields [29] 

£(x, a) = pF p (ai, ■ ■ ■ ,a p ; 61, ■ ■ ■ , b p ; -xa)e a . (15) 

It is worth noting that the weight factors associated with the above CS, play an important role in 
the study of photon number statistics in quantum optics, where similar type of states also appear 
[30]. 
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III. CONNECTION WITH POTENTIALS 



We now proceed to specific potentials for the purpose of illustration and establishing connections 
with the various CS obtained so far. We use the Morse potential [6-9] as an example, for the CHG 
class of potentials, and the Poschl- Teller (PT) potential [11-14] for the HG class. 

The one dimensional Morse potential is given by 

VM = d[l-exp(-ay)} 2 , (16) 

a and d (depth of the potential) being constant parameters. Introducing dimensionless pa- 
rameter = \/2md/ah, dimensionless coordinate q = ay and using the transformation rule: 
x = 2fiexp(—q); the Schrodinger equation yields the eigenfunction 

^{x)^e~ x ' 2 x x l 2 L x n {x), (17) 

where we have set A = 1\x — 2n — 1 for the sake of convenience. Multiplying Eq. (8) from the 
left, by the ground state of the Morse eigenstate: ip^(x) = exp(— x/2)x A//2 , and noting that, for 
b = A + 1 the CHG series can be expressed in terms of the Laguerre polynomials [31] 

$ ( _ n;A+ l ;x) = ___ L A (x) . (lg) 

the coherent state for the Morse potential can be written from the results of the previous section: 

= iv(/3) _1 r(A + i)v M E m + „4.iv (19) 

n=0 l(A + n+ ij 

After normalization one finds, 

Vh{2\x\) 1 (A + n + 1) 
l/3|V2 (/?r A /V e -*/ 2 J A (2v/^). (20) 



Here, I\ is the modified Bessel function of the first kind and J\ is the Bessel function of the first 
kind. The same expression was obtained in [20] using the SUSY based ladder operators. 

In order to see the difficulties associated with properly defining an algebraic structure using 
SUSY based ladder operators [20] 



A^fix) 



x^- ± ^(2n + A + 1 -x) 
dx 2 



(21) 



as alluded to earlier, we explicitly write down the action of the same on the wave functions: 
A + ^(x) = y/(n + l)(n + A + l)ip™ +1 (x) and A~tp^(x) = - ^Jnin + X)tp^(x). In earlier works 
without explicitly giving the diagonal operator, its action was inferred from [A + , A~]%l^(x). How- 
ever, as was noticed in [8] this approach faces problem in constructing the Barut-Girardello CS. 
To better appreciate the difficulties and its resolution, we first identify the corresponding ladder 
operators at the level of monomials through similarity transformations, which keeps the algebraic 
structure intact. These operators act as 

(K+ + D - n)x n = x n+1 (K_ + D-n)x n = n(n + a)x n - 1 . (22) 
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It can be noticed that as compared to our ladder operators, at the level of monomials, the above 
ones contain an additional operator D — n, which yields zero when acting on monomial x n . It can be 
straightforwardly seen that the SUSY based n dependent operators do not lead to a proper algebra, 
a difficulty noticed in [8]. However, the n independent operators lead to the diagonal operator 
K3 = xd/dx + 6/2; together these form a closed SU(1,1) algebra. Hence, for the construction of 
AOCS with a well defined algebraic structure, it is imperative to use n independent operators. 
It can be seen that for the ground state, from which the AOCS are constructed, the above n 
dependent operator is absent. Hence the expression for the CS derived earlier [20] and the one 
obtained here, based on the SU(1,1) algebra, are identical. 

As mentioned earlier the SU(1,1) generators written above are not unique e.g., the following 
three generators also form a SU(1,1) algebra [3]: 



2 d n ,. _,. r d , d 

L + = x - — h 2(A + l)x, L_ = — , and L3 = x- — h 

dx ax dx 



(A + l) 



(23) 



The eigenstate of L_ is given by exp(/3x). The corresponding AOCS, at the level of the wavefunction 
for the Morse potential, can be obtained by making use of the following identities, 



exp 



d 2 /> „s d 



d__ d_ 

dx dx ^ 



d 2 

- X dx-*- {X + 1) Tx 



(24) 



and 



cxp 



x d 2 _ / A + 1 )A 
dx 2 dx 



exp(/3x) = £/3"L^(x). 



(25) 



n=0 



The resulting coherent state turns out to be the Perelomov coherent state for this potential [6]: 

(l-|/3| 2 ) (A+1)/2 „„A/2_.._/ X1 + /3 



^ ER (x,/3) 



yrXATT) 



-x x ^ 2 exp 



21-/3 



(26) 



Hence, the above SU(1,1) algebra gives the operator, whose eigenstate is the Perelomov coherent 
state; recourse has been taken earlier to more complicated nonlinear algebras for this purpose [28]. 

We now derive the CS for the PT class of potentials and concentrate primarily on the sym- 
metric PT (SPT) and PT potentials. Plots of the weight factors associated with the CS of the 
above mentioned potentials will be given along with the quantum carpet structure and the auto- 
correlation figures. The quantum carpet and the autocorrelation plots, transparently bring out the 
phenomenon of revival and fractional revival. 

The SPT potential is 



h 2 a 2 p(p - 1) 



Vspt(v) 



2m cos 2 ay 

and the corresponding eigenvalues and eigenfunctions, in the variable x = sin ay, are 



(27) 



E 



SPT 



2^,2 



h a 
2m 



{n + pf 



n 



0,1,2, • 



a(n\)(n + p)T(p)T(2p) 



_^T{p+l/2)T{n + 2p) 



1/2 



(1 - x 2 y' 2 Ce n {x) 



(28) 



7 



The Gegenbauer polynomials are related to the HG series via the relation [31] 



2 F 1 (-n,n + 2p;p+l/2;z) 



(2p) n 



C£(l - 2z) . 



(29) 



Multiplying (1 — x 2 ) p / 2 from the left in Eq. (10) and using Eqs. (28) and (29), the coherent state 
for the SPT potential is found to be 



Xspt(z,7) = N il) 1 

n=0 

The normalization constant, N^)^ 1 , is 
where, 



Y{2p)T{p + l/2)v^F 



a(n!)(n + p)r(p)r(2p + n) 



1/2 



7>f T (^)- 



(30) 



«r(p) 



s{\i\) = £ 



r( P + i/2)r(2^)0F5(| 7 |) 

| 7 |2n ! 



1/2 



-n!(n + p)r(2p + n) | 7 | 2 p Jq 



2|7l 



dxh p -i{x) 



(31) 



It is worth pointing out that the CS for the SPT potential can be expressed in terms of the Bessel 
functions by using a generating function of the Gegenbauer polynomials [31]: 



Xmpt(x, 7 ) = N{ 1 y 1 Y{p + 1/2)6*° ( \ 



1/2-/9 



^p-l/2(7Vl -x 2 



Similarly we can construct the CS for the PT potential. The PT potential is 



2„,2 



V PT (x) 



h a 
2m 



«(«-!) , pip - 1) 



+ 



sin 2 ax cos 2 ax 



, K,p>l, 



(32) 



(33) 



whose energy eigenvalues and the eigenfunctions are 

h 2 a 2 



2m 



K + p + 2n) 2 , n = 0, 1,2,- •• 



V£ T (x) = C n (cosax) A (sinax) K P^ 1 / 2 '^ 1 / 2 )(l-2sin 2 ax) 
Here C n is the normalization constant and given by [32] 



2a(K + p + 2n)r(n + 1)F(k + p + n) 
r( K + n + l/2)r(p + n+l/2) 



1/2 



(34) 



(35) 



Using Eq. (34) and Eq. (35) in Eq. (10) the AOCS for the Poschl- Teller potential can be written 
in the form 



n=0 



r(K + l/2)(A + l/2) T 



2a(« + 1/2)„(k + p + 2n)r(n + l)r(/s + p + n) 



1/2 



^J{x). (36) 



In the above equation relation between Jacobi polynomials and HG series has been used [31]. As 
has been done for the Morse potential one can also construct Perelomov type coherent state here. 



FIG. 1: The density plot of |^(a;,i)| 2 for, 7 = 10 and the maximum value of n — 20 in the symmetric 
Poschl- Teller potential. Darkness displays a low and brightness a high functional value. 




FIG. 2: Squared modulus |spt (x,"f,t 7^ | x,-f,t — 0)spt | 2 (y-axis) and time (x-axis), for 7 = 10, 
p = 2, and n = 20, of the autocorrelation for the SPT potential. The peaks show the revivals whereas the 
intermediate ones fractional revivals 



The quadratic nature of the eigenspectra of PT and SPT lead to the possibility of revival and 
fractional revival in this quantum system. Keeping this in mind, we now proceed to study the 
time evolution property of the above CS. As expected, these states show a very rich structure 
involving revival and fractional revivals. We give in Fig (1) quantum carpet representing the time 
evolution of the above state. The plots for the auto correlation are provided Figs (3) and (4), 
which clearly bring out the above features. Interestingly, there have been some recent proposals to 
use the fractional revival for the purpose of factorization of numbers [34]. In the above quantum 
carpet the ridges and the valleys follow a curved path, unlike the square-well case where these are 
straight lines [25]. We also notice richer structure arising due to interference. The origin of these 
structures in the square-well case have been understood, the present scenario needs a thorough 
understanding. 

As is well-known, the weight factors associated with the CS carry physical significance e.g., 
these factors for the oscillator CS give rise to a Poisson distribution. The weight factors associated 
with the HG class of CS, derived above, are related to the HG distribution in probability theory 
[33], as is clear from the plot of the weighting distribution in Fig (5). 



IV. CONCLUSIONS 



In conclusion, we have developed a general algebraic procedure for finding the annihilation oper- 
ator coherent states for a wide class of quantum mechanical potentials. Interestingly, the Perelomov 
type CS also emerged, as AOCS, from different realizations of relevant symmetry algebras. Crucial 
use was made of a simple method for solving linear DEs, which gives a precise connection between 
the solution space and the space of monomials. Ladder operators, corresponding to various symme- 
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FIG. 3: Squared modulus |pt (x, 7, t ^ | x, 7, t = 0)pt | 2 (y-axis) and time (x-axis), for 7 = 10, k = 2, 
p = 6, and n = 20, of the autocorrelation for the PT potential. 




FIG. 4: Plot of the weighting distribution |c„| 2 (y-axis) in Eq. (30), with respect to n (x-axis). This plot 
captures the HG statistical nature of the coherent state. Here 7 = 10, p = 2, and the maximum value of n 
is 20. 



try algebras can be identified straightforwardly. The method is potential independent and enables 
one to find the AOCS and connect them with earlier oscillator based approaches. It is applicable 
to quantum problems having infinite number of bound states as well as the ones possessing finite 
numbers. Generalizations to analogs of the nonlinear CS for oscillators is also made transparent. 

Under time evolution, these CS showed revival and fractional revivals. This manifests in the 
quantum carpet structure as well as the auto-correlation functions. Interestingly, these phenomena 
in the context of square- well [25] have led to a proposal for factorizing numbers [34]. The intricate 
structure of quantum carpet needs careful analysis in light of recent proposals to use CS for 
quantum information storage [35]. The weighting distributions associated with these CS also needs 
to be studied more elaborately, in the complete parameter range, for manifestation of non-classical 
behavior. 

Also as a continuation of the present work, it would be interesting to study the features of Wigner 
quasi-probability distributions for these CS, in light of the interesting results obtained recently in 
this area [36]. It is worth noting that, recently the Wigner distribution for the Morse eigenstates 
have been studied [37]. Since the method used here also applies to many-body interacting systems 
it is worth constructing and studying the corresponding CS [38]. A number of these questions are 
currently under study and will be reported elsewhere. 
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